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O . Abstract 

We generalize our previous 2-dimensional model in which a pairing con- 
PLhI densate {ipip) was generated at large A^. In the present case, we allow for 

Oh. both {ipip) and a chiral condensate {'ipip) to exist. We construct the effective 

>-Ch ■ potential to leading order in jj, and derive the gap equations at finite den- 

sity and temperature. We study the zero density and temperature situation 
analytically. We perform the renormalization explicitly and we show that the 
physics is controlled by a parameter related to the relative strengths of the 
interactions in the pairing and chiral channels. We show that although a solu- 
tion to the gap equations exists in which both condensates are non- vanishing, 
the global minimum of the effective potential always occurs for the case when 
one or the other condensate vanishes. 
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I. INTRODUCTION 



In a previous paper [1], we introduced a variant of the Gross-Neveu model which, in the 
large N hmit, exhibits the formation of a pairing condensate (ipip). In that work, we derived 
and solved the gap equation, demonstrated that the coupling was asymptotically free, and 
discussed some of the properties of the relevant Green's functions. Our work was motivated, 
at least in part, by recent interest in the formation of similar condensates in QCD in the 
presence of a chemical potential. [2-4] 

But in real QCD, as well as in a variety of condensed-matter systems, the pairing con- 
densate is not the only possible one. In QCD at low density and temperature there is a 
chiral condensate; it is only as the density increases that one believes the {qq) condensate 
disappears and a new phase, characterized by a non- vanishing {qq), replaces it. Likewise, in 
other systems, which condensates exist depends on external parameters like the density and 
temperature, and also on the relative strengths of various couplings. With this in mind, we 
extend our previous work by considering a more general model governed by two independent 
couphngs: the original Gross-Neveu term [5] that promotes the condensation of and 
the term considered in our earlier paper that produces a (ipip) condensate. As in ref. 1, we 
shall be able to write down the gap equations exactly to leading order in and to renor- 
malize the couplings, thereby rendering the gap equations finite. For the general case of 
non- vanishing density and temperature, we shall be able to write the gap equations in terms 
of a single integral over a momentum variable k; the effective potential is then expressible as 
a further integral over the auxiliary fields. To solve the gap equations will require numerical 
evaluation of these integrals, which takes us beyond the scope of the present work. 

If we set temperature and chemical potential to zero, we can derive a closed- form ex- 
pression for the effective potential, and can explicitly perform the renormalization of the 
couplings. We then find that there is one dimcnsionlcss parameter, independent of the 
renormalization scale, whose value determines which of the condensates is present. This 
situation might be described as "partial dimensional transmutation": the unrenormalized 
theory has two bare couplings whereas the renormalized one has a renormalization scale, 
which is arbitrary, and a dimensionless parameter, independent of this scale that controls 
the physics. We find that the gap equations have three types of solution: two in which one 
or the other of the condensates vanish, and a third, mixed case, in which both condensates 
are non-vanishing. It turns out, however, that the true minimum of the effective potential 
is always at a point where one of the condensates vanishes. Thus the mixed case is never 
realized physically, at least for zero temperature and density. 
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In the remainder of the paper, we shall describe the analysis that leads to the above 
results. In section II, we define the model and derive the gap equations for the general case 
of non-vanishing density and temperature. In section III, we specialize to the case /i — T — 0, 
and derive an explicit form for the effective potential. We renormalize the coupling constants 
and derive thereby the renormalized effective potential and gap equations. Section IV is 
devoted to an analysis of the gap equation and we derive the conditions under which one or 
the other condensate dominates. Section V contains some conclusions. 



II. GENERAL CONSIDERATIONS 



We begin with the Lagrangian: 

+ 2G'^{iij^^j5^^^^){il^^^j5^^^) - (2.1) 

The flavor indices, summed on from 1 to iV, have been explicitly indicated. The first quartic 
term is the usual Gross-Neveu interaction, whereas the second such term, which differs in 
the arrangement of its flavor indices, induces the pairing force to leading order in In the 
flnal term, ji is the chemical potential. 

Strictly speaking, a condensate cannot form, because it breaks the U{1) of Fermion 
number and hence violates Coleman's theorem [6] . Similarly, as well as conden- 
sates cannot exist at flnite temperature in one spatial dimension because of the Mermin- 
Wagner theorem [6]. Nevertheless, it is meaningful to study the formation of such conden- 
sates to leading order in as explained in ref. [7]. 

Our conventions are: 7° = cti; 7^ = — «C2; 75 = (J3. The pairing term, proportional to 
G^, may then be rewritten: 

2G^i^^h^i^^H^h^i^^'^ = -G^[eap^P^^HT][e^si^?i^^P] ■ (2.2) 

Following standard techniques [8] we add the following terms involving auxiliary flelds m, 
St, and B: 

A £ = -^[m + g'i^^Pf - ^(^^ - G^e./3V'i^^V?^^)(S + G\si^^^Hf) ■ (2.3) 

This addition to C will not affect the dynamics. In £' = £ -|- A>C, the terms quartic in 
fermion flelds cancel, and we have 
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™2 of ID 

C = ^{t ij-m- /.7°)^ - ^ - + - ^^o.^i^ ■ (2-4) 

We integrate out and -0^ to obtain the effective action depending on the auxihary fields 
m, B and S^: 

re//(m,5,Bt) = j ^4^(_^_^)_^rrlog^^-'^yl-^Trlog[l + M2(^^)-V2^-V2] (2.5) 

where we have subtracted a constant (independent of the auxihary fields) and have defined: 
^ = 7°(^ i7- - IJ'l^) = ^^0 + ^CTs^s - A* - "i(7i (2.6) 

so that — —ido — icr^dx — A* — mai. 

Since we are looking for a vacuum solution, we have assumed in (2.5) that B,B^ and m 
are constants and have set — AB'^B. The trace on fiavor indices will give a factor N. 
The large- A?" hmit is achieved by setting g'^N — A, and G'^N — k/4, and letting N ^ oo 
with A and k, fixed. We define the effective potential Veff via 



eff 



-N{ I d'x)V,ff (2.7) 



and we therefore have 

V^ffim, M) = ^ + ^ + V}j}{m, M) , (2.8) 



with v}fj:{m,M) = ^[tr\og{A'^A)xx + tr\og{l + M^{A'^)-'^a2A-^a2)xx], where now the trace 
is only over the spinor indices. 

We next generate the local extrema of Vef f by solving 



9Veff ^ dVeff 

dm? dM^ 



. (2.9) 



We evaluate the matrix products in Vg^j in momentum space, with — >■ ik^. The traces 
can be done with the help of 

tr[ ^] = (2.10) 

for any Vq, V. After some manipulation, equations (2.9) become 
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1 _ -^K// _. [ (fk [kl - hi + //2 + M'- 



m 
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^ (2.11 



2 A dm^ J (27r)2 D 

K~ dM^~'J(27ry D ^^■^^> 

where D — [kl — kl — + — — 4{m'^kl + //^/ci — m^kf]. In this expression, ko is 
shorthand for /cq + iesgnko, where e — > 0"^. This prescription correctly implements the role 
of II as the chemical potential. 

The equations can be reduced further by doing the ko integral. Let us define k± — 
V&i ± 262, where bi ^ M'^ + m'^ + /j,'^ + kj, and 62 = [M2m^ + //^(A;^ + m^)]i Then evaluating 
the ko integral by contour methods, taking proper account of the ie prescription mentioned 
above, we find 

— = — dki[— + — + , ^ i^-^)] 2.13 

2A SttJ-a k+ k_ ^MV + iJ?{kl + rn?) K k_ 

and 

K Stt i-A ^ A;+ A;_ ^M'^w? + iJ,^{kj + m^) k+ k^ ' 

The /ci integrals are logarithmically divergent and we have regularized them by imposing a 
cutoff A. This will be absorbed in the renormalization process to be described in the next 
section. Note, however, that the combination jx ~ k given by a convergent integral. This 
fact will ultimately lead to the renormalization-scale independent constant mentioned in the 
introduction. 

We observe from the form of equations (2.11) and (2.12) that the function V^jj- can be 
reconstructed by integrating with respect to and Af^ in the expressions for and ^. 
This will determine Vg^j- up to a single constant V^jj:{0,0), which can be chosen arbitrarily 
without affecting any physical quantity. Explicitly performing this integration we obtain for 
the unrenormalized determinant correction to the effective potential 

V^^'^ (m, M) = -^ dki [k+ + k_] (2.15) 
27r Jo 

To generalize this discussion to the case of non-zero temperature, one returns to eqns. 
(2.11) and (2.12), and one continues to Euclidean space via the replacement ko —ik^ with 
k^ now considered real. The statistical-mechanical partition function is obtained from the 
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Euclidean zero temperature path integral by integrating over a finite regime in imaginary 
time r — it from to P — Because of the cyclic property of the trace, the Fermion 
Green's functions are anti-periodic in r and one has the replacement 

dk,^-Y: (2.16) 

P n 

where the antiperiodicity gives the Matsubara frequencies: 

= k,. = (2.17) 

To do the sum over the Matsubara frequencies, one uses the calculus of residues to obtain 
the identity: 

|E/(iu;„) = -^tanh^ Resfiz,) (2.18) 

where Zg are the poles of f{z^ in z in the complex plane; Resf{zs) is the residue of j{z) at 
Zg and we have assumed the function j{z) falls off at least as fast as l/l^^l^"*"^ for large 
It will be convenient to use: 

tanh^-^ — \ — 2nf{zs) 



z 



where 



n^{z) = 



e^^ + 1 

is the usual Fermi-Dirac distribution function. 

Rotating equations (2.11) and (2.12) into Euclidean space as described above, we get: 

1 dV^]] r dH: [-ki - kl + + - m^] 



2 A dm'^ J (27r)2 D 

r (Pk 
J W) 



1 _ dV}}l _ r (Pk [-kl-kl-i^^-M^ + m^] 
K ~ 9M2 ~ J (27r)2 D 



where now 



where D = [—kj — k\ — + rn? — /i^]^ — 4[—m'^kl + fikf — m'^kj]. There is no longer any 
need for an ie in the definition of k^. Performing the sums over the Matsubara frequencies 
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we obtain the unrenormalized form of the equations which are given by the same expression 
as the zero temperature ones found earher, with the replacements: 

f -^(l-2n/(M) 

f -^(l-2nKM) (2.20) 

As before we can integrate this to get the determinant correction to the effective potential 
which in unrenormalized form is: 

V^^\m,M)^-^ dh[k+ + k^ + l\og{l + 6-^''+) + llog{l + 6-^''-)] (2.21) 

ZTT Jo p p 



III. THE CASE /X = T = 



Renormalization of the effective potential is best discussed in the context of the zero 
temperature and density sector of the theory where we can define the renormalized coupling 
constant in terms of the physical scattering of Fermions at a particular momentum scale. 
This vacuum sector is interesting in its own right and we shall be able, by analytic means, 
to derive the result that depending on the relative strengths of the 2 couplings the theory 
will be in one or another broken phase but never in a mixed phase. Setting fi ~ T — Q we 
obtain 

dV^ff _ 1 r,. M . 1 1 



^ „ . dfci[(l + — ) , +(1 ) , (3.1) 

dm? At: Jo m ^f,2 + + m)^ ^ ^kf + (m - Mf] 



= / dk^lil H ) h (1 ) (3 2) 

dM^ An Jo + (m + M)2 ^ ^A;? + (m - M)^] 



which is solved by 

yW(m,M) = -;^ dki[Jkl + (M + m)2 + Jkf + (M - - 2ki] . (3.3) 

ZTT Jo * 

This can be integrated to give the unrenormalized effective potential: 

1 ?A 2A 
_ [(M + m)Hn{— ) + (M - m)Hn{— -)] . (3.4) 
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We renormalize by demanding that the renormahzed couphngs kr and \r satisfy 



m=mo Kr 

and 



^ - ^ (3.5) 



I M=Mn — 1 — • I""'-*/' 



'^eff 



dw? m=mo Xr 

Here M — Mq, m — mo designates an arbitrary renormahzation point on which the couphngs 
will depend. Using these conditions to solve for A and k in terms of Xr and k,r yields the 
renormahzed form of the effective potential: 

Veff = m^[a + —In \ |] + M^[b + —In \ |] 

47r 7o 47r 70 

^ , M + m , , 

+ —mMln — 3.7 

27r ' M-m ' ^ ' 



where a and h are the following constants: 

1 3 



a — 



2Xr 47r 



6^J--^ + ^^/n|^°-"^°| (3.8) 
kr 27r Stt Mo Mo + mo 

and 7o =| Mq — mg |. 

Note that the renormahzation we have just performed at /i — T — is also sufficient to 
remove all divergences from the effective potential in the more general case of non- vanishing 
chemical potential and temperature. The addition of /i and T will only result in finite 
corrections to the gap equations and therefore to the vacuum values of m and M. We shall 
return to this point in section V. 

The gap equations are properly derived by differentiating Veff with respect to B and m 
and then setting these derivatives to zero. Because Veff depends only on B^B and m^, it 
will always be possible to have solutions with one of m or S or perhaps both set to zero. 
The gap equations can be written 

1 1 I M^ -m^ L M , , M-m , , 

m 2a + — + —In- -] In \ — = (3.9) 

^ 27r 27r 70 ^ 27r ' M + m ' ^ ^ 



and 



1 1 iM^-m^L m, ,M-m, 

M[b + — + —ln^- -]-ir^n — =0. 3.10 

^ 47r 47r 70 ^ 47r ' M + m ' ^ ^ 
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The solutions m — m* and M — M* will give us the local extrema of K//- The first of 
these equations is an identity if m = 0, and the second if M = 0. Also, the values of m 
and M that solve these equations are physical parameters that must be independent of the 
renormalization scale 70. Thus these equations tell us how a and b individually run with 70. 
We note, however, that if we solve for the combination 

5^a-b^ —\— - —]ln I I (3 n) 

the scale 70 drops out. Therefore 5 is a true physical parameter in the theory; we shall see 
in the next section that its value controls which of the two condensates m and M can exist. 



IV. ANALYSIS OF THE GAP EQUATIONS 

It will be useful in the following to note that, at a solution of the gap equations (3.9) 
and (3.10), the effective potential takes the simple form 

K//(m,M) = -i-(m^ + M2) . (4.1) 

Our goal is to analyze all the solutions of the gap equations and to find the one that produces 
the global minimum of K//. This will then represent the true vacuum of the theory. 

There are four types of solution to (3.9) and (3.10). The first is simply to set m = M = 0, 
leading of course to \^ = 0. Clearly, from (4.1) we see that if any other solution exists, V = 
cannot be the minimum of V. The second and third types are obtained by setting M — 0, 
m 7^ and m = 0, M respectively. If M = 0, then from (3.9), we have 

= 7o exp[-{l + 47ra)] (4.2) 

so 

Vo{m, M = 0) = -^e-(^+^'^") (4.3) 
47r 

(we shall use Vq to denote values of V^ff at solutions of the gap eqn.). Likewise, if m = 
0, M 7^ 0, then from (3.10) 

= 70 exp[-{l + ATib)] (4.4) 

Vo{m = 0,M) = e-(^+^"') . (4.5) 
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Thus we see that 



Vo{m = 0, M) < Vo{m, M = 0) if 5 > (4.6) 

and 

Vo(m, M = 0) < Vo(m = 0, M) if 5 < . (4.7) 

The fourth case is when both m and M are non-vanishing. It is then convenient to define 
p = ^ and to combine the gap equations in the form 

^=(/-i)i^+j^+j^M '"'';^'" )] (4.8) 

and 

5= . (4.9) 

An p P + 1 

Both these equations are even in p, so we may take p > for convenience. Eqn. (4.9) tells 
us immediately that if 5 > 0, < p < 1, and if 5 < 0, p > 1. Furthermore, the r.h.s. of 

(4.9) is bounded between and 77-. Hence we conclude: If I 5 |> 77- there is no solution 
with both m and M non- vanishing. If | 5 |< there is such a solution, with the property 
that m > M if 5 > and M < m if (5 < 0. 

It remains to decide whether Vo{m, M) can be the global minimum. To this end, it is 
convenient to reexpress the gap equations once more in the following form: 

- (1 + Aira) = In * ~ * - ln{ I in (4.10) 

7o P + 1 

and 



- (1 + 4nb) - In U _ in{ \ P — 1 |^} . (4.II) 

To P + 1 

Prom these, making use of eqs. (4.1), (4.3) and (4.5), we immediately obtain 

l/o(m = 0, M) = e-(i+^'^'') = 9i{p)Vo{m, M) (4.12) 



and 



yo(m, M = 0) = e-(i+4-'') = ^2(p)l^o(m, M) (4.13) 



where 
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(1 + p) 



p 



(4.14) 



l + p2 



and 



92{p) 



p 



(4.15) 



l + p2 



Eq. (4.12) is the relevant comparison if > S > and < p < 1, whereas eq. (4.13) is 
relevant for > 5 > —77- and p > 1. 

We observe, however, that 5'2(^) = Qiip), so both cases reduce to the following: if we 
can show that gi{p) > 1 in the range < p < 1, then Vo{m, M) is never the global minimum 
(recall that the VqS are all < 0). On the other hand, if gi{p) < 1 in this range, it will be 
possible to have Vo(m, M) be the global minimum. 

To settle this question, write gi — e'*, with 



In the range of interest, p^ < p, so the r.h.s. is a sum of positive terms. Hence h{p) > and 

giip) > 1- 

We conclude that the global minimum of Veff has M = 0,m 7^ if 5 < 0, and m — 
0, M if 5 > 0. 



Prom eqn. (2.21) we can see that the corrections due to non- vanishing temperature and 
density do not affect the ultraviolet behavior of the integrand in the ki integral defining 
V'^^\ Therefore, the renormalization that we have performed at p = T = in section III 
suffices to remove the ultraviolet divergences from the effective potential, and will allow us 
to send the cutoff to infinity. It is perhaps worth recording the complete result explicitly. 
We find, from eqns. (3.5) and (3.6), that 



h{p) = (1 + -) In (1 + p) + (1 - -) In (1 - p) - In (1 + p^) 



P P 



= IYT^^ + ^ (1 + p) + (1 - ^) In (1 - p) . 



(4.16) 



V. CONCLUSIONS 




(5.1) 




(5.2) 
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where a and b are defined by eqn. (3.8), and X is a divergent integral given by 
1 

= T- / dkii- 
Att Jo 



kl + (mo + Mo)2 ^kl + (mo - M^f 
1 2A 

= — [log ( )] + terms which vanish as A ^ oo . (5.3) 

27r VTo 



Thus the full renormalized effective potential may be written 

1 2 

Veff = aiw? + ctaM^ / dki[k+ + A;_ + -log (1 + e~^''+) 

27r Jo 3 



+ ^log + 



- 2A;i - , + , . 5.4 

2 ^fcf + (mo + Mo)2 + (mo - Mo)' 

where ai = ^(1 + 47ra) and 0:2 = ^(1 + 47r6). If cci < a2, then at = T = the vacuum 
has m^ = mp = 7oe^^'^"i g^j^^j ^2 _ g fjgj-g m^ is the dynamically generated fermion 
mass. It is convenient to choose the renormalization scale so that nip = 70. This entails 

a-i = 0,a2 > 0. Furthermore, we are free to choose Mq = 0, so that mo = mp- Then l^// 
takes the form 

1 t^^ 2 
K// = «2M2 - — / dk^[k^ + A;_ + -{log (1 + 6"^*^+) + log (1 + e"^*^")) - 2k^ 
Zn Jo p 

-{m' + M') / ]. (5.5) 

' + nip 



We observe that if we set M = in this expression, we obtain, with E = y/W+m?, 
yeff{m\T,^) = —[log ^ - 1] 

- 1 r s + + ^^-^^ 

which is the effective potential for the Gross-Neveu model in agreement with refs. [9] and 
[10], and furthermore at T = the integral can be done to give an exphcit form that agrees 
with the results of ref. [11]. 

Similarly, in the opposite case 02 < cti, we have, in the fi = T = vacuum, m' — 
and = = 706"^'^"^, where A is the dynamically generated gap. So we choose 
a2 = 0, CKi > 0, and mo = 0, A^ = 70 = Mq . The effective potential becomes 



11 



1 roo 2 

V^ff = aim^ - — / dh\k+ + k^ + -(log (1 + 6"^'=+) + log (1 + e-^''-)) - 2h 

ZTT JO p 
-1 



- (m' + M^) , ^ ] ■ (5.7) 



When = 0, this expression gives us the effective potential at finite temperature for 
the pure Cooper-pairing model considered in ref. 1. Explicitly we have 

^^// = -^^^""^ -^^-pL V^^^ + e-/^^/^] . (5.8) 

Note that it is independent of the chemical potential, as was the case at T = 0. For T » M 
we can expand the integral in eq. (5.8) to obtain [12] 

V^ff = ^[log{^)-^], (5.9) 

where 7 is Euler's constant. The minimum of this function occurs at M = 0, which means 

that the condensate vanishes for large T, as expected. This feature is also borne out by 
numerical evaluation of eq. (5.8) [13]. 

In this paper, we have derived the general forms for the effective potential in leading 
order in eqs. (5.4)-(5.6). We have analyzed the case fi — T — in detail, showing how 
the phase structure is governed by the relative magnitude of the two constants ai and 0:2. 
For II — T — this structure is remarkably symmetric in the two condensates m and M. We 
expect this symmetry to disappear in the general case, however, because n acts to suppress 
m. To study this will require careful numerical analysis of the integrals in eqs. (5.5)and(5.7), 
in order to see how the extrema of Vgff depend on the parameters ai,a2, /i and T. This work 
is currently being actively pursued. It is our hope that the results may be useful for a variety 
of condensed matter systems, for QCD, and perhaps may have cosmological implications as 
well. 
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